﻿ Transformarea z http://shannon etc upt ro/teaching/ps/Cap10 TransfZ pdf Transformarea z generalizeaza transformarea Fourier in timp discret pentru intreg planul complex la fel ca si transformarea Laplace care generalizează transformarea Fourier pentru intreg planul complex Pentru variabila complexa se foloseste notatia: j z x j y r e x y     ; , 22 z r x y   argz  1 intrare Sistem iesire nn 0jnnn h[n]            0 0 0y n h n z z H z 00x n z r e Raspunsul unui sistem discret liniar si invariant in timp la exponentiala complexa discreta de modul neunitar 22j   z x jy re x, y R z r x y arg z ; ; si            k 0jnn n n n k n                  z r e y n h n z h k z z h0 0 0 0 0 ; 0kz kk       kn             H z h k z y n z H z ; 00 k   nn x n c z y n c H z z          k k k k k  kk Pentru orice sistem discret LIT n -functie proprie 0z 2  -valoare proprie  0Hz 1 Transformarea z bilaterala Transformarea z bilaterala a semnalului x[n] discret aperiodic:  nj ; , 0x n z X z x n z z r e r      Z     n Legatura dintre transformata Z si transformata Fourier in timp discret:  n j n n    Z F ;       0x n z x n r e r x n r                   n   j  1z r , x n e x n  Z F                |z|=1 n j n n    X z r x n e r x n            nn      1n r x n l    Transformata z evaluata pe cercul unitar, de raza |z|=1, 3 este transformata Fourier in timp discret Exemple n 1 , 1 x n a n a    n 1nn X z a z az    00nn 11 Daca 1 , az X z z a     1 1az 1 X  j 1ae Regiunea de convergenta (ROC) = valorile variabilei complexe z pentru care X(z) exista X(z)-functie rationala, zerourile = radacini ale numaratorului; polii = radacini ale numitorului Constelatia de poli si zerouri (PZC) 4 2 n 2 1x n a n         11n  n   11nn X z a z az               1        1n n naz n 11    11   0naz a z 11 Daca 1 atunciXz      1  11az 1az  1 az 1 ; a z 1 1az 5 Proprietatile domeniului de convergenta al transformarii z bilaterale 1 Domeniul de convergenta al transformatei bilaterale z nu poate contine nici un pol al acesteia 2 Daca este un semnal cu suportul de durata finita ,x n X z ,   transformata sa bilaterala exista in tot planul, cu exceptia evenual a punctelor 0 si/sau zz   3 Daca este xn un semnal cu intindere infinita spre dreapta si  cercul zr0daca este din domeniul de convergenta al transformatei bilaterale, aztunci toate punctele din plan , cu zr exceptia eventual 0 a punctului de la infinit, sunt din domeniul de convergenta 6 3 4 Daca este xn un semnal cu intindere infinita spre stanga si daca  0cercul ,este inclus in domeniul de convergenta al transformateizr z ntual atunci toate punctele din plan ,0bilaterale cu exceptia evezr a punctului din origine sunt din domeniul de convergenta 5 Daca este xn un semnal cu intindere infinita atat la stanga cat si la  ta al si daca cercul este0dreapta inclus in domeniul de convergenzr transformatei z bilaterale,atunci acesta este o coroana circulara ce include   12x n x n x n 0cercul zr Semnal x[n] cu intindere spre ROC a lui X(z) este 1Semnal x[n] cu intindere 2 dreapta (cauzal); ROC in afara intersectia celor spre stanga (anti-cauzal); cercului de raza R- +7 ROC disc cu raza R doua ROC (cerc) Exemplu nn  3 , 10 , nanxaanx , 1 n  1 , nxnxnxnanx 212 n1  ;  azna 1 az1 n  n1  nna11  a  11 ,  z a11 z1 a 2 zan11 ;  zaa aa1   zaz  a  8 4 Trasformarea z inversa Fie un contur inchis inclus in domeniul de convergenta al transformatei z bilaterale 1kn : se inmulteste cu si se integreaza pe cercul :X z X z x k z z      k 11n n kdz   z dz x k z dz x k  Pentru orice functie olomorfa    k1n kkz   2nfzjd    DC, si orice curba inchisa din DC avem: nnf z dz f z1in , !ndz za  a f z un punct interior curbei Functia 1 este olomorfa in intreg planul complex,  deci si in DC Punctul 0,este interior cercului considerat a 2j     11,1nfznk          1 ; DC0!ndz x n X z z dz 2j   0,zank  z 111n ; Z X z n X z z dz x n x n X z       2j 9 Definirea transformatei z bilaterale prin constelatia de poli si zerouri Daca este o fractie rationala:Xz  M   0kzz 1k ; X z k z DC  N   pkzz 1k suficienta cunoasterea polilor pk kzz si a zerourilor 0Este pentru a cunoaste X z k cu exceptia constantei  Daca se cunoaste in plus si X z zvaloarea intr-un anumit punct  00 din DC, se poate determina si constanta k 10 5 Exemplu Fie constelatia de poli si zerouri din figura 0cu 0 si pzz 0,5 1z      k X z 1: 1Pentru 0,51 0,5zz  Daca semnalul xn este cauzal atunci   DC : zz 0,5 p Deoarece cercul unitar este in interiorul DC semnalu l are si transformata Fourier in timp discret:xn  j 1e  XΩ   jj   0,5 1 0,5ee Fie punctul A de pe cercul unitar, plasat la unghiul  lungimea arcului de cerc in radiani   OA                XX ; arg AP11 cazul general, notand cu si cu 0In polii si zerourilepk kzz si cu si unghiurile pe care le formeaza vectorii Az k k pk 0si cu axa orizontala, se poate scrie:kAz M 0kAz MN 1k X z k ;  arg k         Nkk 11kk pkAz 1k Frecventa  reprezinta lungimea arcului de cerc masurat in radiani, pe cercul unitar, de la intersectia acestuia cu semiaxa 0 si pana in punctul , sensul fiind celxA trigonometric 12 6 Transformata z unilaterala  n Z x n z X z x n z    uu 0n Transformata zz unilaterala poate fi privita ca si transformata bilaterala a unui semnal cauzal Domeniul de convergenta al transformatei unilaterale este fie tot planul, fie exteriorul unui disc cu centrul in origine Transformata z unilaterala este indicata pentru studiul sistemelor cauzale caracterizate de ecuatii cu diferente finite liniare si cu coeficienti constanti care au conditii initiale nenule 13 Proprietatile transformarii z Notatii: uZZ x n X z z DC x n X z   , ;   xu Tema de curs: uZZ Demonstratiile , ; y n Y z z DC y n Y z     yu 1 Liniaritate ; , cel putinax n by n aX z bY z z DC DC      xy     uuax n by n aX z bY z 14 7 2 Translatia in timp 0n       0 , x n n z X z z DC Demonstratie    0n n m     00mnnnm                Z x n n x n n z z z x m z00 n m m          00Daca 0, 0 se elimina din DC iar daca 0 , se elimina din DC n z n z In cazul transformarii unilaterale, pentru 0 :n 0  0n n m  0nnm             00uZ x n n x n n z x m z z   00n m n 1   0nmm ,z x m z x m z         00m m n  1   0nn           0ux n n z X z x n z 15 0nn  3 Modularea in timp   00j n j , e x n X e z z DC   Demonstratie n     0 0 0j n j n jn Z e x n e x n z x n e z,    nn    00j n j e x n X e z   u Mai general:    nnz z z        00 DC; uz x n X z x n X 0 0 0z z z    4 Reflectarea semnalului 11 , x n X z DC     z Demonstratie m  11mn     n Z x n x n z x m X            nmzz      16 8 5 Diferentierea in domeniul n 1  1 1 ; ;x n x n z X z z DC            1 1 1 1 x n x n z X z x               u  6 Insumarea in domeniul n nXz   z **  ; DC DC ; DC 1 x k X z z z z           1  11kzz   1     uX z x k n k  ; 1 x k z    1  1kz  Demonstratie m 1 Fie y n x 1 1 ;k x n y n y n X z z Y z                        k  1 Y z X z     1 1z Pentru transformata unilaterala: 1 117 1 1 1 , 1 x n y n y n X z z Y z y y x k                            u  k  7 Diferentierea in domeniul z dX z  ; DC nx n z z    dz  udX z nx n z  dz 8 Conjugarea complexa in domeniul timp * * * , DC x n X z z   * * * x n X z   u Demonstratie * n   * * * * *n ;Z x n x n z x n z X z        nn     18 9 9 Teorema convolutiei in domeniul timp , x n y n C      DC DC cel putin x n y n X z Y z z          xy x n y n X z Y z         uu Demonstratie  n Z x n y n x n y n z                n   nk    k x k y n k z z      nk      n km  mk  y m z x k z Y z X z          mk      19 10 Teorema produsului semnalelor in domeniul timp 1z du , DC x n y n X u Y       2j u u      DC ; DCz R z R z R z R        x x x y y y z     si R u R R R     x x y y u Prin inmultire           DCR R z R R z R R z R R            x y x y x y x y  1z du , DC x n y n X u Y     u  2uj u u  20 10 Demonstratie  11n n n    Z x n y n x n y n z X u u du y n z                    2j  nn     n  11z du z du      X u y n X u Y             22j u u j u u      n    Particularizari * 21z du 2 * *                x n l y n x n x n X u X *   2j u u  2 * 1n  *z du x n z      * X u X   2nj u u     *j j j    Daca si siu e u e du je d     2 j   *1jjed         jx n X e X e j   n-2je  22  11j  21 X e d X d          22     11 Teorema valorii initiale aunui semnal discret cauzal 0 lim lim x X z X z   u zz  Demonstratie 1n  0 1 X z X z x n z x x       u 0nz Trecand la limita se obtine relatia din enunt 12 Teorema valorii finale a unui semnal discret cauzal lim 1 lim 1 x z X z z X z       11uzz Demonstratie   1nm 1mnn   1 1x n x n z x n z X z x m z X z             uu 0 0 1n n m     m 0 1 0 z x m z x X z z X z zx         uu 0m   n Pentru 1: lim 1 limz x n x n z z    1 0 ,X z x     u 11zz 0n sau: 0 lim 1 0 lim 1 x x z X z x x z X z           22 11uuzz 11 Relatia dintre transformata z si transformata Laplace Prin esantionare se obtinex t X s  aa  ˆ       a s sx t x nT t nT n  snT ˆ sx t x nT t nT x nT e    LL     a e s a e nn   nn        a s d d d a sx nT x n Z x n x n z x nT z nn  ; sT sTx t t Z x n x t t Z x n   L L    ssa T d a T u du  sse z e z Tr Laplace a semnalului analogic esantionat Tr Z a semnalului digital corespunzator x t t Z x n x n x nT  ;      a a sTddsT ss ez sT sze23 Studiul sistemelor discrete liniare si invariante in timp prin intermediul transformarii z Forma simpla a teoremei convolutiei semnalelor discrete face din transformata Z un instrument util pentru studiul sistemelor discrete LIT Daca sistemul studiat este cauzal si are conditii initiale nenule, este utila transformarea Z unilaterala 24 12 Functia de sistem a unui sistem discret, liniar si invariant in timp Functia H(z), numita "functia (de) sistem" sau "functia de transfer a sistemului" caracterizeaza complet comportarea sa în planul complex z, dupa cum raspunsul la impuls h[n] caracterizeaza complet comportarea sistemelor in timp 25 Daca sistemul discret este stabil, exista si transformata Fourier in timp discret a lui h[n] si prin urmare cercul unitar este in domeniul de convergenta al functiei de sistem ce caracterizeaza un sistem stabil Daca sistemul este cauzal: H(z) = H(z), atunci u domeniul de convergenta al functiei este exteriorul unui disc Daca sistemul este stabil si cauzal, cercul unitar este in DC Toti polii functiei sistem a unui sistem stabil si cauzal au modulul subunitar ceea ce este ehivalent cu faptul ca sunt cuprinsi în interiorul discului unitar: |z| <1 p26 13 Determinarea raspunsului unui sistem discret liniar si invariant în timp, utilizand transformarea Z Daca x[n] este dat si se specifica h[n] ↔H(z) atunci se aplica semnalului de intrare transformarea directa rezultand X(z) Raspunsul in complex este Y(z)=H(z)X(z) Aplicand transformarea Z inversa, rezulta semnalul raspuns, y[n] Daca se lucreaza cu sisteme si semnale cauzale în conditii initiale nenule, se va utiliza numai transformarea Z unilaterala 27 Calculul transformarii z inverse Sunt aplicabile trei metode: 1 Calculul direct al integralei, 2 Transformarea functiei Y(z) intr-o suma de fractii simple, 3 Metoda dezvoltarii functiei Y(z) in serie de puteri 28 14 2 Transformarea functiei Y(z) intr-o suma de fractii simple Metoda se aplica in cazurile in care Y(z) este o functie rationala Functia Y(z) poate fi un -1 raport de polinoame în z sau z -1 Recomandam sa se lucreze în z, notand -1 pentru comoditate z = x, deoarece in majoritatea cazurilor tabelele sunt date în -1 functie de puterile lui z 29  xRxN 1  xIzYxY  xDxD kk   knczcxcxI kkk kkk  baxRks kim  i 1    xxxxxDmki km  xR   xxa;  mxxmm  xD  kis  xRds1  k  xxb    kkxxkiski  xDdx!is   k Se aplica transformarea inversa fiecarui termen din suma, dupa ce se revine mai intai la z-1 = x, utilizand tabelele de transformate si tabelele de proprietati ale transformarii 30 15 ks a m kibRx   i   1mikmD x x xxx  k si Rx    1ksRxd k   mma x x     sikkib x x Dx   !kDxsidx    k  mxx kxx 31 Exemple z  1 zzzYzY50DC u  z z25050 1 xz16888  zY; 11  xxxxzz424242 44168  zY 1111 z,z,zz250150142 nn    n,,ny250504  2 transformAceeasiata ,zzzY250DCdar nn    n,,ny1250504  3 transformAceeasiata ,z,zzY50250DCdar nn   n,n,ny25041504 32 16 3 Metoda dezvoltarii functiei Y (z) in serie de puteri Dezvoltarea functiei Y(z) in serie de puteri in jurul originii Exemplul #a 1z , : 0Y z e ROC z    1 1 21 1 1 1nnz     1e z z z z        01! 2! ! !nnn 1 yn   !n 33 Exemplul #b z , : 0,Y z e ROC z z       21 1 1 1mmz 1e z z z z        01! 2! ! !mmm 01 11  nn 1zz    !!nn     nn       11       1y n n n n            n !1!nn      34 17 Exemplul #c 1 ln 1 , :Y z az ROC z a    semnal cauzal     11nnnn  11aa     n     ; 1Y z z y n n       1nnn Teorema valorii initiale: 1 0 lim ln 1 0y az       z n 1   n 1y n a n       n 35 Exemplul #d 1n Y z z a y n a n,           1 1az 36 18 Exemplul #e, aceeasi transformata ca si la #d, ROC diferita 1z Y z z a  ,  1 1zaaz semnal anticauzal, transformata z contine doar puteri ale lui z 1 1mn n n 1Y z a z a z y n a n               1mn   37 Sisteme discrete liniare si invariante in timp, caracterizate prin ecuatii cu diferente finite liniare si cu coeficienti constanti NM         0; 0kka y n k b x n k a 00kk M k NMkbz  0kkkNz functia de transfera z Y z b z X z H z     kkN  00kkkDz kaz 0k Functia de transfer a unui sistem discret LIT: functie -1 rationala in z or z Zerouri : radacinile numaratorului N(z); Poli : radacinile numitorului D(z) 38 19 Pentru un sistem cauzal si stabil: polii sunt in interiorul cercului unitar |z|=1 1z pk Gr M Teorema valorii initiale: Nz  0 lim lim lim - finitah H z H z      u z z z   Dz  Gr N Gradul numaratorului functiei de transfer a unui sistem cauzal si stabil este mai mic sau egal cu gradul numitorului, M N 39 Contributia polilor unui sistem discret cauzal in raspunsul la impuls al acestuia 40 20 2 Pereche de poli dubli complex conjugati,   *ppjj z r e z r e si p p p p Contributia lor la functia de transfer este de forma: ** 1 1 2 2a a a a        2211 ppjj    11ppjj 11r e z r e z 11ppr e z r e z     pp Contributia lor in raspunsul la impuls este de forma: nn               12sin sin p p p p p pA r n A nr n n  Daca r 1, raspunsul este amortizat-sistem stabil p Daca 1, raspunsul creste in timp-sistem instabil r p 41 Calculul raspunsului unui SLIT discret caracterizat printr-o ecuatie cu diferente finite Se considera SLIT cauzal, cu conditii initiale nenule (dar semnalul de la iesire este necauzal): NM      , 0 0kka y n k b x n k a 00kk NM a Z y n k b Z x n k     k u k u 00kk NkMk  k n k n        k u k ua z Y z y n z b z X z x n z  01k n n01k  Daca semnalul de intrare este cauzal, 0 pentru 0:x n n    N k M  k n k a z Y z y n z b z X z        k u k u  0 1 0k n k    Cunoscand si conditiile initiale, se determina si apoi X z Y z y n42   uu 21 Exemplu nj 0                   1 , , ,ynkenxnxnayny01 k  1             1 ; zzyzYazzY1 juu1 0 ze1   ayk1   zY   ju111   0       azazze111 j 0    aykake111      jjj111   000     azazeazeea111   njn11  0 keka n 1               1 ; 1 , zanyany1  jj  00  eaea   43 Sisteme de ordinul I kzk              zzHnkxnayny 1 ; DC  1 1  azaz OA 1jj          n           h n ka n H e e PA PA Daca 0 1a maximul raspunsului in frecventa pentru minim, 0, filtru trece jos PA   Daca 1 0,a   maximul rezulta pentru , filtru trece sus  44 22 Sisteme de ordinul doi kzk2   21 zHnkxnyanyany 21  azazzaza1221 2121 aaa42 Avem211 zeroun de ordinul 2 originein si polii z ,p21 2  Daca2 4 cei doi poli sunt conjugaticomplex ezaaj p21 ,21 Daca2 4aa polii reali sunt 21 Se considera sistemul cauzal si se cauta conditiile pe trebuiecare sa le scaindeplinea iicoeficient si aaea stabilitat apentru obtine 21 aaa4222 1121 ; 4aaa 212 2 a2 1 , aa:realisunt1 dacasau polii 22 4 24222 ; aaaaa-04 21211 Rezolvand sistemul de iinegalitat :rezulta 042 ; 1 ; aaaaaa1 121221 45 Parabola corespunde existentei unui polsingur real dublu  j21  Raspunsul frecventain al :sistemului ekH21 APAP 21 46 23 Functia de sistem echivalenta unor sisteme discrete conectate in serie si in paralel       nhnhnh;       nhnhnh; e21 e21       zHzHzH       zHzHzH e21 e21 47 Forme de implementare ale filtrelor numerice utilizand transformarea Z NMMN  , ,1 knyaknxbnyNMaknxbknya  0kkkk 1000kkkk Forma specifica de reprezentare a implementarii in 48 forma directa I S-a considerat a = 1 0 24 Forma directa I (N+M sumatoare, fiecare cu 2 intrari) Forma specifica de reprezentare a implementarii in forma directa I Un sumator cu 2N intrari S-a considerat M=N, a=1 49 0 Forma directa II (2*N sumatoare, Forma specifica de reprezentare a fiecare cu 2 intrari) S-a considerat implementarii in forma directa II M=N Doua sumatoare fiecare cu N intrari S-a considerat M=N, a0=1 50 25 